Polymer Welding: Strength Through Entanglements 
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Large-scale simulations of thermal welding of polymers are performed to investigate the rise 
of mechanical strength at the polymer-polymer interface with the welding time t w . The welding 
process is in the core of integrating polymeric elements into devices as well as in thermal induced 
healing of polymers; processes that require development of interfacial strength equal to that of the 
bulk. Our simulations show that the interfacial strength saturates at the bulk shear strength much 
before polymers diffuse by their radius of gyration. Along with the strength increase, the dominant 
failure mode changes from chain pullout at the interface to chain scission as in the bulk. Formation 
of sufficient entanglements across the interface, which we track using a Primitive Path Analysis is 
required to arrest catastrophic chain pullout at the interface. The bulk response is not fully recovered 
until the density of entanglements at the interface reaches the bulk value. Moreover, the increase 
of interfacial strength before saturation is proportional to the number of interfacial entanglements 
between chains from opposite sides. 



Thermal welding is a common means of joining poly- 
meric elements (TJ [5] . Two polymer surfaces are brought 
into close contact above their glass transition tempera- 
ture T g and are allowed to interdifTuse for a welding time 
i w . Polymer melt dynamics suggests that for homopoly- 
mer melts the properties of the weld should be indis- 
tinguishable from those of the bulk once the chains have 
diffused by about their radius of gyration i? g . In practice, 
however, welds reach the bulk strength at much shorter 
times [2J. A key question is what determines the rapid 
rise in interfacial strength and how it is related to mass 
uptake across the surface and polymer entanglement. En- 
tanglements are topological constraints of polymers by 
other chains [3] that control their bulk visco-elastic and 
plastic response. 

Experiments have quantified the strength of welds by 
measuring the interfacial fracture toughness in tensile 
fracture and the peak shear strength of lap joints 0]- 
6 . Both quantities grow as tl/ 4 at short times and then 
saturate at their bulk values. Several different molecular 
mechanisms have been proposed to explain the scaling of 
strength with t w . Some assume the strength is simply 
related to interpenetration depth [7J-|n| or the areal den- 
sity of chain segments bridging the interface [JO] or the 
contour length of bridging segments [11] . In many cases 
these models are motivated by physical pictures of the 
development of entanglements at the interface, but chain 
friction may also be important |12j . 

It is difficult to distinguish between the proposed 
strengthening mechanisms with experiments. Entangle- 
ments are not directly observable and experiments are 
usually restricted to post-analysis of fracture surfaces [T3] 
or bulk scattering [8] , which does not isolate the response 
of the failing region. In contrast, computer simulations 
provide full spatial resolution throughout the failure pro- 
cess, allowing macroscopic stresses to be directly related 
to molecular structure and dynamics. Recently devel- 
oped methods also enable tracking of entanglements on 



a microscopic level [I4HT6] . 

In this Letter, we present results from large-scale 
molecular dynamics (MD) simulations of welding be- 
tween surfaces of highly entangled homopolymers. The 
interfacial strength after a welding time i w is determined 
by a simple shear test that mimics experiments [3H2]- 
As in the experiments, the interfacial strength rises lin- 
early with tl/ 4 before saturating at the bulk value well 
before the time for the chains to diffuse a distance 
_R g . The dominant failure mechanism changes from chain 
pullout at small f w to chain scission at large f w and in 
the bulk. Evolution of entanglements during welding is 
tracked using a Primitive Path Analysis (PPA) algorithm 
[Til ITT] . The crossover to bulk response coincides with 
a rise in the entanglement density at the interface to the 
bulk density. Moreover, the areal density of entangle- 
ments between chains from opposing surfaces is linearly 
related to the interfacial strength at small t w . 

Our simulations employ a canonical bead-spring model 
|18j that captures the properties of linear homopolymers. 
Each polymer chain contains N spherical beads of mass 
Tn. All beads interact via the truncated shifted Lennard- 
Jones potential 

U L j(r)=^u [(a/r) 12 -(a/rf-(a/r c ) 12 + (a/r c ) 6 } , (1) 

where r c is the cutoff radius and Uu(r) = for r > r c . 
All quantities are expressed in terms of the molecular 
diameter a, the binding energy uq, and the characteristic 
time t = c^m/ito) 1 / 2 . 

For equilibration, beads along the chain were con- 
nected by an additional unbreakable finitely extensible 
nonlinear elastic (FENE) potential 

f/ F ENE(r) = -^ J R ln[l-(r/ J R ) 2 ] , (2) 

with i?o = l-5a and k = 30woa~ 2 . For mechanical tests, 
chain scission plays an essential role and a simple quartic 
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potential was used 

U Q (r) = K(r-R c ) 2 (r-R c ){r-R c -B) + U , (3) 

with K = 235lM /fcs, B = -0.7425a, R c = 1.5a, and 
Uq = 92.74467ito- This potential gives the same equi- 
librium bond length as C/fene and prevents chains from 
crossing each other so that entanglements can be studied. 
The bonds break at a force that is 100 times higher than 
the interchain C/lj, consistent with experiments and pre- 
vious simulations [19[I20], Previous work has shown that 
the entanglement length for this model is N e = 85±7 and 
that the mechanical response for N = 500 is characteris- 
tic of highly entangled (large N) polymers |2"TH2"4"] . 

The equations of motion were integrated using a 
velocity- Verlet algorithm with a time step St < O.Olr. 
The temperature was held constant by a Langevin ther- 
mostat with a damping constant T |18| . A million r will 
be abbreviated as IMt. All simulations were carried out 
using the LAMMPS parallel MD code [2"S] . 

Two thin films were constructed following the stan- 
dard methodology discussed by Auhl et al. [2B]. Each 
film contains M = 4800 chains of length N = 500 beads 
or a total 2.4 million beads. Periodic boundary condi- 
tions were applied along the x- and y- directions with di- 
mensions L x — 700a and L y — 40a. The thickness in the 
^-direction was maintained at L z = 100a using two repul- 
sive confining walls. Each film was well equilibrated at a 
temperature T = LOuo/fce with r c = 2.5a, T = O.lr -1 
and pressure P = maintained by expansion/contraction 
along the x-direction. To form the welding interface at 
z = 0, the films were placed as close to contact as pos- 
sible without overlap. The interdiffusion was along the 
z-direction. During interdiffusion, volume was held fixed 
by repulsive walls perpendicular to the z-direction. 

After welding for a time i w , the system was quenched 
rapidly below the glass temperature T g sa 0.35uo//cs. 
First the cutoff radius was reduced to r c = 1.5a to de- 
crease computational cost, reduce density changes and fa- 
cilitate comparison with past mechanical studies [HO |2"1T 
[24] . Then the temperature was quenched at constant vol- 
ume with a rate T = — 1Q^ 3 uq/ {Ubt) to T = 0.5moAb 
where P = 0. Subsequent quenching to T = 0.2iio/fcs 
was done at t = -2 x 10 _4 -ito/(fcsr) and P = 0. A Nose- 
Hoover barostat with time constant 50r was applied to 
P xx and P yy . The repulsive walls were maintained at 
z = ±L Z . We verified that our conclusions are not sensi- 
tive to the details of the quench protocol or geometry. 

Shear was applied to the interface in a manner similar 
to a shear test of a lap joint [2J. Beads within 5a of the 
top and bottom were held rigid and displaced at constant 
velocity in opposite directions along the y— axis. The 
average strain rate in the film, drf/dt = 2 x 10 _4 t _1 , 
was low enough that it did not affect the mode of failure 
and stress had time to equilibrate across the system [27] . 
The shear stress a was determined from the mean lateral 
force per unit area applied by the top and bottom walls. 
The temperature was maintained at T = 0.2uo/fcB with 



a Langevin thermostat (T = It 1 ) acting only on the 
x-component to avoid biasing the flow. 
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FIG. 1. (a) Shear stress a versus shear strain 7 for the 
indicated t w . (b) Average interpenetration depth (d) versus 
strain, (c) Number density pbroken of broken bonds versus 7. 
The inset shows the maximum shear stress <r max normalized 
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by its average bulk value <r m " x versus t 



Stress-strain curves for polymer welds at different t w 
are compared to the bulk response in Fig. [ija). The 
bulk curve is typical of amorphous polymers. A narrow 
elastic region is followed by yield and strain hardening - 
a gradual increase in stress with strain. At sufficiently 
large strain, the material begins to fail. The shear stress 
reaches a maximum value cr^^ and then drops. 

Even for the shortest t w , the stress-strain curve fol- 
lows the bulk response up to 7 ~ 4. As t w increases, the 
stress follows the bulk curve to larger 7. For t w > 2Mr 
the response is nearly indistinguishable from bulk. In 
experiments, the strength of the weld is characterized by 
the maximum shear stress cr max before failure and plotted 

against tl/ 4 to test scaling predictions [2]. Our results, 
shown in the inset of Fig. |T](c), are very similar to ex- 
periments [2] . There is a constant strength from van der 
Waals (Lennard- Jones) interactions at short times before 
appreciable interdiffusion. The strength saturates to the 
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bulk value as t w increases past lMr and is statistically in- 
distinguishable from bulk behavior by about 3Mr. This 
is much shorter than the time for polymers to diffuse by 
their radius of gyration id 3> lOMr [55] ■ At intermediate 
times there is a linear rise in strength that is consistent 
with cr max ~ tl/ 4 . Both simulations and experiments are 
limited to a factor of ~ 3 change in stress that prevents 
precise testing of power law scaling. However the tl/ 4, 
scaling is motivated by an assumption that strength is 
proportional to the interdiffusion distance. This is di- 
rectly tested as described below. 

Simulations allow the change in maximum stress to be 
correlated with changes in failure mechanism and molec- 
ular conformation. Strain leads to stretching and orien- 
tation of polymer chains that is directly related to strain 
hardening [23] |29l [30]. Starting near 7 ~ 4 we see a 
noticeable tension along polymer backbones that grows 
rapidly with 7. This tension acts to pull back any chain 
segments that have diffused across the interface. One 
way to quantify chain pullout is to measure the evolu- 
tion of the average interpenetration depth (d) of beads 
that have crossed the initial interface (z = 0) as a func- 
tion of strain (Fig. [TJb)). For t w < O.lMr, (d) begins to 
drop at the strain where the stress deviates from the bulk 
response and decreases most rapidly near <7 max - The final 
interfacial width (d) is only the size of single bead and 
corresponds to complete chain pullout. At this point the 
stress reaches a constant value that represents the friction 
between two separated films. 

At larger t w , chains cannot be pulled out from the op- 
posing surface and (d) decreases by at most a single bead 
diameter. The chains have diffused far enough that the 
tension required for chain pullout is high enough to break 
bonds along the backbone. Fig.[ljc) shows the strain de- 
pendence of the number density Pbrokcn of broken bonds 
averaged over the entire volume between moving layers. 
The curves for t w > 2Mr are similar to the bulk. Bonds 
begin to break above 7 = 6 and the highest rate of bond 
breaking is reached near 7 = 10. This point coincides 
with the peak in shear stress that indicates mechanical 
failure. For O.lMr < t w < 2Mt there is a crossover where 
the number of broken bonds rises rapidly and the amount 
of chain pullout measured by the drop in (d) decreases. 
At smaller t w , bond breaking is localized near the initial 
interface. For t w > 2Mr bonds break uniformly through- 
out the system, confirming that the interface has become 
as strong as the bulk. It is interesting to note that a sim- 
ilar transition from chain pullout to scission occurs with 
increasing chain length in previous simulations of craze 
formation JTHl HU [22] . Chain scission only occurs when 
chains are long enough to form entanglements (typically 
N > 2N e ) that prevent chain pullout. By analogy, it is 
natural to expect that the transition to scission and bulk 
response occurs when chains have interdiffused enough 
to form entanglements at the interface. 

Entanglements have proved elusive in experimental 
studies. However, the representation of entanglements as 
binary contacts between the primitive paths of polymer 
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FIG. 2. Density profiles for (a) total and (b) interfacial TCs 
at the indicated i w - The inset shows the reduced maximum 
shear stress versus the areal density of interfacial TCs N^ c /A. 
Error bars are comparable to symbol size. 



chains has provided many insights into the properties of 
polymer melts |14H16j . The primitive paths are obtained 
by fixing the chain ends and minimizing the chain length 
without allowing chain crossing. To limit excluded vol- 
ume effects, the chain diameter is then decreased by a 
factor of 4 [T7] . Contacts between the resulting primitive 
paths are counted to determine the number of topological 
constraints (TCs) . We found that the ratio of the density 
of TCs, ptCj to the bulk density, /o!^ k , was insensitive to 
the precise details of the procedure used to identify TCs. 
Past bulk studies show that pxc is proportional to the 
entanglement density [TH - TTT] and we refer to TCs and 
entanglements interchangeably below. 

Figure [2ja) shows the normalized density of TCs as 
a function of height relative to the initial interface. At 
small f w the chains have not interdiffused enough to pro- 
duce any entanglements at the interface, but there are 
two peaks near z = ±10a. These peaks reflect the fact 
that polymers near free surfaces are compressed perpen- 
dicular to the surface [3TJ [35] . Chains in this pancake- 
like anisotropic conformation are subject to more TCs. 
As welding proceeds, diffusion increases the density of 
entanglements at the interface and reduces the peaks on 
either side. By 2.5Mr the density has become nearly 
uniform across the system. 

To correlate between cr max and entanglements more 
quantitatively, we focus on the number of TCs between 
chains that are on opposite sides of the interface at 
t w — 0. All of these formed by interdiffusion and it is nat- 
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ural that these interfacial entanglements should be most 
important in strengthening the interface by preventing 
chain pullout. Fig. [2jR) shows the profiles of the normal- 
ized interfacial TC density p^c- ^ s ^ ne polymers diffuse, 
interfacial TCs spread outward from the interface and 
grow in number. The insert in Fig. [2] (b) shows the corre- 
lation between the normalized weld strength ^max/^max 
and the areal density of interfacial TCs, N^ c /A. There 
is a linear correlation between strength and interfacial 
entanglements at short times. At long times the inter- 
facial strength saturates while the number of interfacial 
entanglements continues to grow. 
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FIG. 3. The interfacial density of TCs, N\ c /A versus the 
average interpenetration depth (d) and a linear fit. 

From the chain-packing model of entanglements (331 — 
135] , the number of interfacial entanglements N\ c scales 
with the volume spanned by chains that have diffused 
across the interface. This can be estimated as 2A(d). 
Figure 3 confirms that there is a linear relation between 
N\ c /A and (d) over the entire range studied [36]. As 
noted above, welding models have generally assumed that 
&max rose linearly with (d) and then used reptation the- 
ory to argue that both scale as tl/ 4 . Subsequent work has 
shown that interfacial diffusion is more complicated be- 
cause of the anisotropy that produces the peaks in entan- 
glement density noted above and the prevalence of chain 
ends at the interface Indeed, recent simulations 

show that there is not a single simple scaling exponent 
over the range of times studied here [28 . Figures 2 and 



3 show that whatever the time dependence is, the funda- 
mental factor determining strength is the entanglement 
density. One of the main differences between welding 
models is that some have assumed a simple proportion- 
ality between N\ c /A and (d) EH E2] , while others 
assume a minimum interpenetration distance is needed 
for entanglements [9] [37] [3B] . Fig. 3 supports the latter 
interpretation, with a minimum distance ~ 1.5a that is 
about a third of the distance needed for saturation. 

To summarize, we have demonstrated that the devel- 
opment of interfacial strength during welding is closely 
related to the formation of entanglements across the in- 
terface. The interface becomes mechanical indistinguish- 
able from surrounding regions when the bulk entangle- 
ment density is recovered at the interface. There are then 
sufficient entanglements to prevent chain pullout at the 
interface and the joint fails through the bulk mechanism 
of chain scission. Before the bulk strength is recovered, 
the interfacial strength rises linearly with the areal den- 
sity of interfacial entanglements. This quantity is not 
accessible to experiments, but is linearly related to the 
interdiffusion distance which has been measured. These 
findings should help further development of theoretical 
descriptions of entanglement evolution across a polymer- 
polymer interface and constitutive molecular modeling of 
fracture in polymers. Of particular interest will be stud- 
ies that vary the entanglement density by changing the 
entanglement length or making systems immiscible. 
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